Quasi-phase-matching of high-order-harmonic generation using multimode 

polarization beating 
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The generalization of polarization beating quasi-phase matching (PBQPM) and of multi-mode 
quasi-phase matching (MMQPM) for the generation of high-order harmonics is explored, and a 
novel method for achieving polarization beating is proposed. If two (and in principle more) modes 
of a waveguide are excited, modulation of the intensity, phase, and/or polarization of the guided 
radiation will be achieved; by appropriately matching the period of this modulation to the coherence 
length, quasi-phase-matching of high harmonic radiation generated by the guided wave can occur. 
We show that it is possible to achieve efficiencies with multi-mode quasi-phase matching greater 
than the ideal square wave modulation. We present a Fourier treatment of QPM and use this to 
show that phase modulation, rather than amplitude modulation, plays the dominant role in the case 
of MMQPM. The experimental parameters and optimal conditions for this scheme are explored. 
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I. INTRODUCTION 

High harmonic generation (HHG) is a nonlinear pro- 
cess that enables the production of tunable, coherent soft 
x-rays with applications in time- resolved science (H-Q) 
ultrafast holography |4|, or diffractive imaging [5]. An 
important feature of HHG is that it is simple to achieve: 
focusing driving laser radiation to an intensity of order 
10 14 W/cm 2 in a gaseous target yields coherent radia- 
tion with frequencies corresponding to the odd harmon- 
ics of the driving radiation uj. A semi-classical theory of 
this phenomena has been developed by Corkum [(| and a 
uantum treatment has been given by Lewenstein et al. 



However, without additional techniques, HHG is 
highly inefficient - with a typical conversion ratio of 10 -6 
for generating photons of energy in the 100 eV range, 
and 10~ 15 for generating 1 keV photons. This low effi- 
ciency is partially due to fact that the driving field and 
the harmonic field have different phase velocities. As 
a consequence a phase difference develops between the 
driving field and harmonics generated at each point in 
the generating medium; this in turn causes the intensity 
of the generated harmonics to oscillate with distance be- 
tween zero and some maximum value along the direction 
of propagation, z. The phase velocity difference is char- 
acterized by the wave vector mismatch, Afc, which arises 
from neutral gas, plasma, and waveguide dispersion; it is 
given by Afc = k{qu) — qk{u) where q is the harmonic and 
k(uj) is the propagation constant for radiation of angular 
frequency u). The distance it takes for the two fields to 
slip in phase by 7r is the coherence length L c = n/ Afc. 
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One way of avoiding the phase-mismatch problem is 
to balance the dispersion so that Afc = 0, a situation 
we will describe as "true phase-matching" in order to 
distinguish it from the quasi-phase-matching discussed 
below. With true phase- matching - assuming absorp- 
tion can be neglected - the intensity of the harmonics 
grows quadratically with the propagation distance z 0- 
10l ] . With long wavelength drivers for phase matching, 
conversion efficiency can be achieved up to 10~ 3 for the 
VUV region and 10 -6 in the x-ray region [10]. However, 
without long wavelength drivers, true phase-matching 
can only be achieved for relatively low-order harmon- 
ics — corresponding to low photon energies — since at 
the higher driving intensities required to generate high- 
order harmonics the dispersion becomes dominated by 
the free electrons and cannot be balanced by the other 
terms Q. For higher-order harmonics, so-called quasi- 
phase-matching (QPM) may be employed in which the 
harmonic generation is suppressed in the out of phase 
zones, enabling monotonic growth of harmonic intensity 
as a function of z. Techniques for QPM include counter- 
propagating pulses [liTlilil ]. multi-mode beating (l6l - [l9| . 
and modulated waveguides (20| . 

Multi-mode QPM (MMQPM) relies on coupling in two 
or more waveguide modes |16l4la ]. If the two modes 
travel at different phase velocities, then the intensity will 
beat along the propagation length, thereby modulating 
harmonic generation resulting in QPM . In this paper, 
we investigate the effect on HHG of the modulation in 
both the intensity and phase of the beating driving ra- 
diation. We show that under certain conditions, phase 
modulation due to mode beating enables harmonics to be 
generated with greater efficiency than ideal square-wave 
QPM modulation. Moreover we show that MMQPM is 
dominated by modulation of the phase, rather than the 
intensity, of the driving radiation — an effect which was 
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not considered in our earlier analysis (l6l - [l8j . 

Recently we have proposed a new class of QPM based 
on modulation of the polarization state of the driving 
radiation within a hollow core waveguide (2l| - [23 |. Here 
we discuss one example of polarization-control QPM: po- 
larization beating QPM (PBQPM) [HI, [H- In this ap- 
proach, a linear birefringent system modulates the polar- 
ization of the driving pulse between linear and elliptical. 
Because harmonic generation is suppressed for elliptically 
polarized light, QPM can be can be achieved if the pe- 
riod of the polarization beating is suitably matched to 
the coherence length. This paper describes a generaliza- 
tion of MMQPM and PBQPM which combines these two 
schemes: multi-mode polarization beating quasi-phase 
matching (MMPBQPM), which utilizes beating between 
two waveguide modes to modulate the intensity, phase, 
and/or polarization of the guided radiation. These mod- 
ulations can lead to QPM if the coherence length L c of 
the harmonic generation is appropriately matched to the 
beat length, L&, which is the distance it takes for the two 
modes to develop a phase difference of tt. In addition 
to controlling the relative input polarizations of the two 
modes, the relative polarization angle between the two 
modes can also be controlled. This increased parameter 
space affords greater opportunities for QPM. 

In addition we further analyze MMPBQPM using a 
Jacobi-Anger and Fourier decomposition which affords 
additional insight into the processes leading to QPM. 
Similar Fourier techniques used to analyze quasi-phase 
matching for HHG can be found in [25|, |26[ . 

The paper is organized as follows: in Sec. lU the equa- 
tions for MMPBQPM are derived for two modes, Sec 
IIIII presents the the Jacobi-Anger and Fourier analyti- 
cal analysis, and Sec. IIVI discusses numerical simulation 
results. 



Here, the polarizations of the modes m — 1 and m = 2 
are respectively taken to be parallel to, and at an angle 
O, to the x-axis: 



E 1 (r) 



M 1 (r) 




2{T > ~ 1 M 2 (r) sin(fi) 



(2) 



(3) 



where M m (r) is the normalized transverse electric field 
profile for the mth mode. 

At time t' = t + At, where At = —/3i/u)z, the electric 
field components are given by, 
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ciMi(r) cos(W) 




(4) 



( c 2 M 2 (r) cos(fi) cos(wt' + A/3z) 
V c 2 M 2 (r) sin(O) cos(wt' + A/3z) 



where Af3 — f3i — /3 2 , and the beat length Lb = 7r/(A/3). 

From this, the values of t' corresponding to the max- 
imum and minimum electric field amplitudes can be 
found, from which the ellipticity is given by, 



e{r,z) 



l n{€[t' max (r,z);r,z]y 
m[t' mm (r,z);r,z}Y 



(5) 



where 5ft[£] 2 = $l[£ x } 2 
major axis is given by, 



^[e;^] 2 , and the angle of the 



0(r, z) = tan' 



-i ( K[<£vW ma x(r,z);r,z] 
\$t[<£ x ][t' max (r,z);r,z] 



(6) 



II. DERIVATION OF THE ENVELOPE 
EQUATION 

A. Mode Propagation Equations 

In this section, we develop the set of general mode 
propagation equations for two linearly polarized modes 
with azimuthal symmetry. If we assume that two modes 
are excited, then within the waveguide the electric field 
may written in cylindrical coordinates as, 



i(r, z) = E ic 1 E 1 (r)S^ +ia ^ z - Ui ^ 
+ c 2 E 2 {r)e l[{fi2+ ' la2)z - ult] ^ 



(1) 



where r is the radial coordinate from the propagation 
axes, E\ and E 2 are the normalized transverse electric 
fields of the driving mode m = 1 and modifying mode 
m = 2 respectively, E n is the electric field amplitude 
constant, j3 m and a m are the propagation constant and 
damping rate of mode m, and oj is the angular frequency 
of the driving radiation. 



It is useful to note that the relative relative intensity 
of the driving wave at any given point is given as: 

J(r, z) = <E x (r, z)£* x (r, z) + <E y (r, z)£* y (r, z) (7) 

+2c 1 M 1 c 2 M 2 cos(n)cos(A J 8z) (8) 

assuming that the damping terms are small. For the 
remaining of the paper, we will focus on r — and will 
define Mi = Mi(0), M 2 = M 2 (0) to simplify notation. 
Moreover, to simplify arguments, we will stipulate the 
following normalization condition: c 2 M 2 + c|M| = 1 
where we have taken c m M m to be real, which can always 
be achieved by a suitable shift of the z or t coordinates. 



B. Derivation of the growth equation 

If we write the electric field of the qth harmonic as 
1 



E 



Aqui) i[k{qui)z-quit] 



c.c 



(9) 
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then, within the slowly-varying envelope approximation, 
the equation for the growth of the amplitude of the qth 
harmonic becomes, 



2ik(qcu)e lk{qu]) 



dz 



NL 



(10) 



where Pm? is the component of the non-linear polar- 
ization oscillating with angular frequency qui. Now, 
= F'{I,e)e iq ^ i - z \ where F'(I,e) gives the depen- 
dence of the nonlinear response on the intensity and el- 
lipticity of the driving field, and <f> p (z) is the phase of the 
driving field of the pth polarization component. We may 
write 4> x (z) = k\{ijS)z+ij) x {z) and(j)y(z) — k-z{u))z-\-il) v (z), 
where ip p (z) is the additional phase arising from inter- 
ference of the waveguide modes and k m ((d) — f3 m (ui) is 
the waveguide propagation constant for the driving pulse. 
Henceforth all equations will refer to the qth harmonic, 
and so in order to avoid clutter we will drop the (qui) 
superscripts. The growth equation for the amplitude of 
the harmonic for the x- and y- components may then be 
written in the form, 



dz 



9^ =F B (J,e)e- <Afcl V*-W =«p* 



-!Afc 2 Zg#,(z) _ 



F 



3 +i* x (z) 



(11) 

(12) 



where *B X and ^f y are the total phase for the x- and y- 
components respectively; ^} is a projection term that re- 
lates the nonlinear polarization to x- and y- components 
of the envelope function, as developed below; Ak m = 
k(qui) -qk m (u); and F(I, e) = ifi {qui) 2 F'(I,e)/2k(qui). 
We note that, as discussed below, F(I,e) is in general 
complex since the phase of the nonlinear polarization 
depends on the trajectory of the ionized electron, and 
therefore on both the intensity and ellipticity of the driv- 
ing field. In the equation above, we have factored these 
phase terms into ^ x and ^> y . 

Considering first the rr-polarization, the driving field 
may be written as, 

E x (f,z,t) = [dMx + c 2 M 2 cos fte~ lA/32 ] e *fa*-<*) 

= u x ( Z y^ lz -^ (13) 
where A/3 = /3i — fa. Hence we find, 
ip x (z) = arg[w x (z)] 
= tan -1 
c 2 M 2 



-C2M2 cos 51 sin( A/3z) 



ciAfi 



ci Mi + C2M2 cos SI cos(A/3z) 
cos ft sin( A(3z) 



(14) 



where in the last step, the approximation is valid if 
ci ^> I C2 1 ■ Similar considerations show that ip y = 
and (j)y(z) = fa. 

The strength of the nonlinear polarization, e)|, 
depends on the intensity and polarization of the driving 
laser field. Evaluation of \F(I,e)\ requires a model of 



the interaction of the driving field with the atom, as, for 
example, developed by Lewenstein et al. Q ■ However, for 
our purposes it is sufficient to assume that the amplitude 
of the nonlinear polarization, F 
form: 



can be written in the 



F[I(z),e(z)] = AG[I{z)]H{e{z)\ 



(15) 



where A cx ^{qui) 2 /2k(qui) is a constant, G[I(z)} is the 
intensity-dependent term ranging from [0, 1], and H[e(z)] 
is the ellipticity-dependent term ranging from [0, 1]. 

For the purposes of illustrating the operation of MMP- 
BQPM it is sufficient to assume that the intensity- 
dependent term takes the form of a power law G(J) Rj 
i^/ 2 . We will assume rj rj 6, in accordance with earlier 
work [l8|, H3, [H| ; but note that the the broad conclusions 
of the present paper do not depend strongly on the value 
of rj. 

It is also well known that the single- atom efficiency 
of HHG depends sensitively on the polarization of the 
driving laser field [5^] which arises from the fact that 
the ionized electron must return to the parent ion in or- 
der to emit a harmonic photon. Following the argument 
given in [211] , for a given driving intensity the number of 
harmonic photons generated as a function of ellipticity 
maybe approximated by: 



h{e) 



1 



1 



(16) 



where H(e) — yjh(e). 

It is predicted that u = q — 1 within the perturba- 
tive regime, as verified [29J by Budil et al. for harmonics 
q = 11 to 19, and by Dietrich et al. for harmonics up to 
q w 31 [13] ■ Schulze et al. found that for higher-order 
harmonics the sensitivity of harmonic generation to the 
ellipticity of the driving radiation is lower than predicted 
by Eqn (|16|) with fj, = q — 1 [3l|, although in this non- 
perturbative regime the efficiency of harmonic generation 
still decreases with e. Further measurements of the de- 
pendence of harmonic generation on e have been provided 
by Sola et al. [32j |. It is recognized that Eqn (fT6|) is an 
approximation, but it will serve our purpose of demon- 
strating the operation of MMPBQPM. 

The offset angle and ellipticity of the harmonics gener- 
ated by elliptically-polarized radiation have been shown 
to depend on the ellipticity and intensity of the driving 
radiation, and on the harmonic order [23, [H, [U |H, [34| . 
Propagation effects can also play an important role. 
Since the amplitude with which harmonics are generated 
decreases strongly with increasing ellipticity, we are most 
interested in the ellipticity of the harmonics generated for 
small e. It has been shown that for higher-order harmon- 
ics, and/or high driving intensities, both the ellipticity 
and change in ellipse orientation of the harmoni cs g ener- 
ated by radiation with e Rj are close to zero [27]]. We 
will therefore make the simplification that the generated 
harmonics are linearly polarized along the major axis of 
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the driving radiation and that we resolve separately the 
harmonics polarized along the fast and slow axes of the 
waveguide. Thus, the projection term may be written as: 



cos[6(z)] 
sin[9(z)] 



(17) 



and by following the arguments of [21], the coherence 
lengths for harmonics polarized parallel to the x- and y- 
axes are different, and hence for a given Lf, it is only possi- 
ble to quasi-phase-match one of these components. Thus, 
for the remainder of this paper, we will focus on analyz- 
ing the harmonics polarized along the x-axis. Therefore, 
we can approximate I q = t; x (,* + £ y £* s» £ x £*. 

Moreover, the phase of the nonlinear polarization de- 
pends on the intensity of the driving radiation [3^, |36| 
and its ellipticity [13, [34[. We will ignore the effect of 
ellipticity on the phase of F(I, e) since, as shown be- 
low, harmonic generation is dominated by those regions 
in which the driving radiation is close to linear polariza- 
tion. We may write the intensity-dependent phase as a 
Taylor expansion around Iq, 



*(J) = $(Io)+f \ Io (I-Io) + 
« $o + vq (I - I ) . 



(18) 
(19) 



where v = d$>/dI\ lQ . For simulations in this paper, we 
assume v/q w 0.2 rad/10 14 W cm -2 based on previous 
calculations [36[. 

We may now gather the contributions to the total 
phase ^ x {z): 

V x (z) = -$(/„) - Ak lZ + q^ x (z) -qv{I- I ) (20) 

Co Mo 

*o - A/fciz - q— — -cosfisin(A/3z) (21) 
c\M\ 



-1qvc\M\CiMi cos fl cos (A/3z) 



(22) 



where \&o = — ^C^o) — IqvcxMxc-iMi cos fi, and the ap- 
proximation holds if | C2 1 <C |ci|. From Eqns. (fT2")) and 
(j22f , we can rewrite the the differential equation for the 
x-component as: 



^ = C os[e(z)}\F[I(z),s(z 



(23) 
(24) 



III. ANALYSIS OF THE GROWTH EQUATION 

A. Phase Analysis using the Jacobi-Anger 
Expansion 

The exponential term in Eqn (|24[) can be expanded into 
the products of two infinite sums using the Jacobi-Anger 
Expansion: 



e -i*i(z) = e -iAk lZ 



E il M-p> 



+UAfiz 



E M-7> 



VijAfiz 



(29) 



where Ji and Jj are Bessel functions of the first kind 
and l,j G Z. It is insightful to factor terms of constant 
a = I + j to give: 



where 



kiz e ^ eKT 



l+j=tr 



A0z 



(30) 



(31) 



We see that the modulation caused by intensity depen- 
dent phase and mode beating can be resolved into har- 
monics oA/3 of the difference in spatial frequency A/3 of 
the two modes. 



B. Source Amplitude Spatial Fourier Analysis 

The analysis above suggests that it would be useful to 
write the source modulus T as a superposition of Fourier 
components with frequency kA/3 (with k £ Z). For pe- 
riodic modulation of the driving radiation, the source 
modulus can be written as a Fourier scries: 



r(*)= E V ^ K 



A0z 



(32) 



K — — OO 



and hence the growth differential equation can be written 
as: 



where 



1 d£ x e _ lAklZ ^ Utr y Ke Ha+ K )Apz t (33) 



A' dz 



A'{z) = Ae+™° (25) 

T(z)=cos[Q(z)}\F[I(z),e(z)}\ (26) 

= -* x (z) + * (27) 

= Akiz + 7sin(A/3z) + pcos(A/3z) (28) 



in which 7 = g ^Mi cos ^ i s the mode interference 
phase term and p — 2qvc\M\CiMi cosJl is the intensity- 
dependent phase term. 



The terms that contribute to monotonic harmonic 
growth are those for which the phase is stationary, in 
other words, ^f- = 0. This implies that for QPM we 
require Aki = (a + k)A/3. We see that the harmonics of 
the modulation frequency allow QPM of larger wave vec- 
tor mismatch Afc or, equivalently, of shorter coherence 
lengths L c ^\. The fundamental modulation spatial fre- 
quency A/3 has a period 1L\, = 27r/A/3, and hence we may 
write the QPM condition as L\, = (a + k)L c .i = nL Ci i, 
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where n = a + n is the order of the QPM process. Fac- 
toring all the terms contributing to monotonic harmonic 
growth, and ignoring the oscillating terms, the growth 
equation becomes: 



1 d& 
A> dz 



a-\-K=n 



(34) 



for a fixed n, keeping in mind that each of the terms of the 
sum is complex and may have different signs. Eqn (|34j) 
can easily be solved, from which the harmonic intensity 
is found to be: 



1 



A'Sz 2 



(35) 



where S = J2a+ K =n U(tV k . It is useful to note that the a 
and k terms result from phase and intensity modulation 
of the driver respectively. 



IV. SIMULATION RESULTS 
A. Detailed simulations for Lb = 2L c .i and Lb = L c .\ 

To test these ideas, we have conducted a series of sim- 
ulations. Fig. Q] presents the results of simulations for 
Lb = 2L c .i (n — 2) and three values of VI and C2M2 
for q = 51 while Fig |5] presents the same parameters for 
Lb = L C) i (n — 1), ft — 0° and two different values of 
C2M2 for q = 51. These values are compared against ideal 
QPM, which is defined by the square wave modulation 
between zero and one of the harmonic generation with a 
period of 2nL Ct \. 

When O = 0° MMPBQPM is equivalent to "pure 
MMQPM" since the driving radiation remains linearly 
polarized at all points within the waveguide; this is seen 
in FigQJCol (1) where Fig|T]-(l)(b) indicates that modu- 
lation of the source term arises from mode-beating alone. 

When fl = 90°, MMPBQPM is equivalent to PBQPM 
since the mode beating causes the polarization of the 
driving radiation to beat in an analogous way to PBQPM 
driven by a linearly polarized beam propagating in a bire- 
fringent waveguide. This is seen in Fig. [TjCol (3). More 
specifically, as seen in Fig. [TJ- (3)(b), the modulation of 
the source term is seen to arise from solely polarization 
beating. It should be noted that for f2 = 90° the simu- 
lations presented here agree with earlier calculations of 
PBQPM |2l|. For intermediate values of ft (such as in 
Col (2) where fl = 60°), modulation of both the inten- 
sity and polarization of the driving radiation play a role 
in QPM. 

Fig. Q] also compares the growth of the calculated har- 
monic intensity with the approximation of Eqn. (|35p . It 
can be seen that the approximation agrees closely with 
the exact calculation, indicating clearly the dominant 
role played by the terms for which a + n = n = 2a,s 
seen in Fig. [T]- Row (a). 

Moreover, Figfl] - Row (d) maps the values of U a V K , 
modulus phase (the terms in the sum in Eqn (|34p). as a 
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FIG. 1. Simulation results for Columns - (1) "Pure MMQPM" 
(c 2 M 2 = 0.002, n = 0°) (2) "Hybrid MMPBQPM" (c 2 M 2 = 
0.01, n = 60°), and (3) "Pure PBQPM" (c 2 M 2 = 0.05, fi = 
90°); all for Lb = 2L C and q = 51. Row (a) shows the relative 
HHG intensity. Dot dashed cyan line shows ideal square wave 
QPM, solid red line shows the indicated form of MMPBQPM, 
solid black line shows the intensity for no phase matching, and 
the dashed blue line shows the shows the harmonic intensity 
calculated from Eqn (|35[) . Row (b) shows the modulus of 
the intensity-dependent source term G[I(z)] (solid red line), 
ellipticity-dependent source term H[e(z)] (solid green line), 
and the combined source terms \F\ — GH (dashed black line). 
Row (c) shows ^ x {z), the total phase of d^ x /dz (solid black 
line) and ^ x for large z (dotted light red line). Row (d) shows 
the U a V K distribution as a function of a and k for n — a + n 
for n — 2 



function of a and k where a + n = n for a fixed n = 2. 
Hence, only where a + n = 2 is U a V K is nonzero. For the 
case of Pure MMQPM, Col (1), and Hybrid MMPBQPM, 
Col (2), the dominant contributing term is (a, k) — (2, 0) 
indicating that QPM arises predominantly from phase 
modulation of the driver, not intensity modulation. This 
can also be seen in Fig. [TJlc and Fig. [T]-2c where the 



6 



(1) First 
Ressoance 



(2) Second 
Ressoance 



& 15 



O 
K 
EC 




20 



15 



10 




1 

0.5 




FIG. 2. Simulation results for Columns - (1) the first reso- 
nance (C2M2 = 0.0009, Q, = 0°) and (2) the second resonance 
(c 2 M 2 = 0.0091, n = 0°); all for L b = L c and q = 51. Row 
(a) shows the relative HHG intensity. Dot dashed cyan line 
shows ideal square wave QPM, solid red line shows the indi- 
cated form of MMPBQPM, solid black line shows the inten- 
sity for no phase matching, and the dashed blue line shows 
the shows the harmonic intensity calculated from Eqn (|35[1 . 
Row (b) shows the modulus of the intensity-dependent source 
term G[/(z)] (solid red line), ellipticity-dependent source term 
H[e(z)\ (solid green line), and the combined source terms 
\F\ — GH (dashed black line). Row (c) shows ^(z), the 
total phase of d^ x /dz (solid black line) and ty x for large z 
(dotted light red line). Row (d) shows the U a V K distribution 
as a function of a and k for n — a + k for n = 1. 



regions of harmonic growth occur for points where the 
phase ^> x is within 7r/2 of the phase of £ x for large z. In 
contrast, for the case of Pure PBQPM, Fig. Q]- Col(3), 
the dominant term is (a, k) = (0,2). This suggests, as 
expected, that for Pure PBQPM, phase modulation does 
not contribute to QPM but only the modulation of the 
amplitude of the source term caused by caused by polar- 
ization beating. 




FIG. 3. Relative HHG amplitude |£e(;z)| for q = 51 at large 
z (z — 10L Ci i) as a function of Lb/L c ,i and m = 2 mode 



coupling strength c 2 M 2 where cf Mi = 1 — 
to ideal square wave QPM for: (a) £1 = 0°, 



c 2 M 2 , normalized 
and (b) fi = 90° . 



Fig. [5] presents the same parameters in Fig. [T] for 
n = 1 (or Li, = L c \), for two different values of C2-M2 
and = 0°. We see that for both columns, the only 
terms which contribute are those for which a + k = 1, 
as expected. For Col (1), optimal MMQPM enables har- 
monics to be generated with intensities greater than for 
ideal square wave QPM. As indicated in Fig. |2J-(l)(c), 
the region of harmonic growth coincides with ty x being 
within ±7r/2 of the phase of £ x for large z . Moreover, 
the largest contributing term of U a V K in Fig. [2]-(l)(d) is 
(a, k) — (1,0); this suggests that QPM is caused primar- 
ily by phase modulation, and not by amplitude modula- 
tion as reported earlier for MMQPM |Ha - [l8l |. Moreover, 
the phase modulation explains why higher growth than 
ideal square-wave QPM occurs. Fig. [2]- Col (2) shows 
the output at a different mode mix where c^M-x = 0.0091 
and C\M\ — 0.9909. As discussed below, the mode mix- 
tures for which results are shown in Fig [3] correspond 
to two of the peaks in a plot of the output of harmonic 
q — 51 as a function of ci-Mi and C2-M2. 



B. Parameter Space Scans 

This section presents a series of parameter space scans 
for optimizing the harmonic generation by MMPBQPM. 
In an HHG experiment, pressure, coupling angle f2, and 
the mode mix ratio of c\ to C2 are parameters that can be 
adjusted. Pressure tuning equates to tuning the coher- 
ence length, or tuning the ratio Lb/L Cl i assuming that 
Lt is fixed for a specific pair of driving and modifying 
modes. 

Fig. Hi shows, for the MMQPM case (CI = 0°), the 
variation of the harmonic output as a function of Lb/L C} i 
and mode mix C2-M2. Note that here the magnitude of 
the harmonic amplitude, not intensity, is plotted in or- 
der to show more clearly the variation of the harmonic 
output. As expected, MMQPM is optimized for integer 
n. Moreover, the peaks shift to increasing C2-M2 with in- 
creasing n. When n = 1, the QPM condition becomes 
(t+k = l+j + n = 1. The three lowest-order solutions sat- 
isfying this condition are {/, j, k} = {1, 0, 0}, {0, 1, 0}, and 
{0,0,1}. We therefore expect peaks in the HHG inten- 
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FIG. 4. Relative HHG intensity at z — Lb for Lb = 2L Ci i as a 
function of Qmax = tan" 1 ^fjff^ ) and g, normalized to ideal 
QPM. 
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FIG. 5. Relative HHG amplitude for q = 51 after as 
a function of coupling angle O and m = 2 mode coupling 
strength (%M% where c^Mf = 1 — c^Mf , normalized to ideal 
QPM for: (a) L b = L c ,i, and (b) L b = 2L c ,i. 



sity to occur for values of C2-M2 corresponding to peaks 
in |Ji(-p)J (-7)|, |J (-/3)Ji(-7)|, or | J (-p) J (-7)l- 
The maxima along the line Lj,/ L Ct \ shown in Fig[3£i arise 
from the variations of ciMi which optimize the functions 
of \Ji(-p)Jj(-7)\. 

When n = 2, <j + k = l+ j + K = 2, and the three lowest 
order terms are {1,1,0}, {0, 1, 1}, or {1, 0, 1}. Thus the 
values of optimal for C2-M2 will be around the extrema of 

\Ji(-p)Ji(-i)\, |Jo(-p)Ji(-7)l, |Ji(-p)Jo(-7)l- Be- 
cause the positions of the local extremas of the Bessel 
function increase with increasing \l\ or optimal p and 
7 and thus optimal C2-M2 will increase as well. There- 
fore, increasing n will result in larger values of | Z | and 
\j\ contributing to the harmonic growth corresponding 
to the Bessel function peaks shifted to higher values of 
(— p) and (—7) and hence higher values of C2-M2. 

Similarly, Fig. [3}j shows the PBQPM case where 
tt = 90°. As discussed in [Hj], PBQPM will occur when 
Li, = nL c ^i and n is even - as is evident in Fig. [3}d. Since 
fl = 90°, p = 7 = 0, and since J ; (0)Jj(0) = unless 
I = j = 0, monotonic harmonic growth can only occur for 
a = I + k = 0. Hence optimal PBQPM occurs when the 
Fourier coefficient Vk is large for even k. Furthermore, 
PBQPM does not contribute to any phase modulation 



because a = as seen from Eqn (|30|) . The optimal value 
of C2M2 increases with the order n of QPM since increas- 
ing this parameter shifts the Fourier spectrum of the driv- 
ing intensity modulations to higher orders. The optimal 
value of C2M2 is explored more clearly in Fig. |4l which 
shows the normalized harmonic intensity for il = 90° as 
a function of the harmonic order q and the maximum an- 
gle Qmax = tan" 1 ( c^Mi ) ^ ne major axis of the elliptical 
driving radiation makes with the x-axis. If Q max is too 
close to 0°, then the ellipticity modulation is not enough 
to suppress the destructive zones. If Q ma x is too close 
to 90°, then the harmonic generation suppression zone is 
too large to create efficient harmonics. 

Fig. [5] plots, for the cases Lb = L c ^ and Lb — 2L Cj i, 
the calculated relative amplitude at z = 2Lb of the q = 51 
harmonic as a function of f2 and the relative intensity of 
the m = 2 mode. In the case of Fig. Lb = L c .i, the 
relative amplitude achieved with "pure MMQPM" (i.e. 
f2 = 0) is greater than that of ideal square wave QPM 
as explained above. For = 0°, the intensity oscillates 
with increasing C2-M2, with the size of the resonant peaks 
decreasing with increasing C2-M2. These resonance peaks 
are caused by peaks of the products J;(— p) Jj(— 7) with p 
and 7 being a linear function of C2A/2 (for small C2M2) as 
explained above. For the case of pure PBQPM, i.e. ft = 
90° and Lb = L Ci i, the harmonic intensity is seen to be 
very low and almost independent of the relative intensity 
of the two modes since the phase-matching condition for 
lowest-order PBQPM is not satisfied, and QPM is not 
achieved. 

Lowest-order QPM occurs for n — 2, as shown in Fig 
[SJd, a region of bright harmonic generation occurs for 
tt = 90° and c 2 M 2 « 0.04. When ft = 0° the harmonic 
intensity oscillates in a similar manner to that observed 
in Fig. [5^,. 



V. CONCLUSION 

We have developed a generalized analysis of MMQPM 
and PBQPM together with a simplified Fourier analysis 
which gives additional insights into the dominant con- 
tributions of quasi phase-matching. In addition we have 
shown that PBQPM could be achieved without a bire- 
fringent waveguide by exciting a pair of waveguide modes 
with two orthogonal polarizations. 

Our analysis of MMQPM showed, in contrast to our 
earlier analysis 0, [l8| , that QPM is dominated by the 
modulation of phase of the harmonic source term, not 
of its amplitude. This allows, under optimal conditions, 
MMQPM to generate harmonics with an intensity greater 
than possible with ideal, square-wave QPM. 
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